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I. INTRODUCTION 


The causal set program proposes one of a number of approaches to the problem of quan¬ 


tum gravity 1 - [j| . The key aspect of this approach is that it postulates that spacetime is 
discrete rather than continuous. A causal set is a set C of elements ay G C, and an order 
relation such that the set C = {ay, -<} obeys properties which make it a good discrete 
counterpart for continuum spacetime. These properties are that (a) the set is transitive: 
Xi -< Xj -< Xk =>- Xi -< Xk] (b) it is noncircular, ay -< Xj and Xj -< Xi =>• Xi = xf, (c) it is 
locally finite in the sense that the number of elements between any two related elements 
Xi -< Xj is finite, i.e., [ay, x 3 \ < oo; and (d) it is reflexive, ay -< a, V x G C. The action of 
the order relation is to mimic the causal ordering of events in macroscopic spacetime. Since 
all events in spacetime are not causally related, then not all pairs of elements in the set are 
ordered by the order relation. Hence, a causal set is a partially ordered set. 

The continuing development of the causal set program can be roughly divided into two 
categories: kinematics and dynamics. Development of a causal set dynamics is the ultimate 
goal of this program and would represent its version of a theory of quantum gravity. However, 
the ability to investigate the implications of this dynamics will require development of the 
mathematical and computational tools needed to describe the physics of causal sets. This 
collection of techniques make up causal set kinematics. Of particular interest here is the 
correspondence between spacetime as a causal set and macroscopic spacetime. If a causal 
set comprises the true structure of spacetime it must produce a four-dimensional Lorentzian 
manifold in macroscopic limits. 

The goal of determining how we can recognize and extract manifold-like properties from 
a partially ordered set is a kinematic issue of significant current interest. This paper reports 
on a numerical study concerning one aspect of this goal; specifically, we address the cor¬ 
respondence between paths in a causal set and geodesics in continuum spacetimes. To the 
best of our knowledge, this is the first published investigation of this correspondence that 
includes curved spacetimes. In what follows we first describe the key ideas related to this 
correspondence, we then describe the details of how our calculations were performed, and 
finally we present our results and conclusions. 
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II. THE PATH-GEODESIC CORRESPONDENCE 


The properties of causal sets stated previously do not guarantee that the sets will display 
the manifold-like behavior required by the physical theory of causal sets we ultimately seek. 
In fact it is well known that most sets satisfying those properties will not be consistent with 
points in a spacetime manifold. A necessary (but not sufficient) requirement for a causal set 
to be like a manifold is that it can be embedded into a manifold uniformly with respect to 
the metric. An embedding of a causal set is a mapping of the set onto points in a Lorentzian 
manifold in such a way that the lightcone structure of the manifold preserves the ordering 
of the set. The most probable embeddings will be uniform if the mapping corresponds to 
selecting points in the manifold via a Poisson process. When this latter requirement is met, 
the embedding is said to be faithful. 

Recall that the length of the geodesic between two causally related events gives the 
longest proper time between those events. To see what the most natural analog to geodesic 
length is for causal sets we must first define a few terms. A link, A, in a causal set is an 
irreducible relation; so, X; A x*, iff $ Xj 3 x* -< Xj -< Xk- A chain in a causal set is a set 
of elements for which each pair is related; for example, x a -< x b -< • • • -< x 2 _i -< x 2 is a 
chain from x a to x,. A maximal chain, or path , is a chain consisting only of links, such as 
x a ^x b <---< x 2 _]_-d x 2 . 

Myrheim argued >4| that the length of the longest path between two related elements in a 
causal set is the most natural analog for the geodesic length between two causally connected 
events in spacetime. The length of a path is defined to be the number of links in that 
maximal chain. A more formal statement, that we refer to as the Myrheim length conjecture 
is the following: 

Let {C n } be a sequence of causal sets C n = {x*, -<} of increasing density p that 
are faithfully embeddable into a Lorentzian manifold Ad by a map g : C n —> Ad. 

Then, in the limit p —> oo, the expected length of the longest maximal chains 
between ordered pairs (xj,Xj) G C n is directly proportional to the geodesic 
lengths between their images g(xi,xf) G Ad. 

Brightwell and coworkers have proven the above statement for the case when the 
Lorentzian manifold is Minkowski space JfJ. Their result uses the fact that, in ('/-dimensional 
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Minkowski space, the volume V of a region of spacetime bounded by the lightcones of two 
causally connected events is V — mdT d , where r is the longest proper time between those 
events and rn ( ] is a constant. In curved spacetimes (with nonconstant curvature) that re¬ 
lationship breaks down in that (a) any proportionality factor between V and r d will be a 
function of spacetime and (b) the dependence of V on r may include additional powers of r. 
While it is widely believed that the Myrheim length conjecture is true for curved spacetimes, 
no proof has been achieved. 

Underlying the development of the causal set program is a much wider conjecture, termed 
the Hauptvermutung, which effectively says that in the appropriate limits, a physically inter¬ 
esting causal set will bring forth all of the structure of an approximately unique spacetime 
manifold ij ,7']. This conjecture might be stated the following way: 


Let there be a random process that produces a faithfully embeddable causal set 
C. Then, if 3 maps ()\ and g 2 3 g\ : C —■► Adi and g 2 : C —> _Ad 2 , then M.\ and 
A4 2 are approximately isometric. 


One only requires an approximate isometry to allow for random fluctuations where “approx¬ 
imately isometric” might mean that the manifolds are close in the sense discussed in [ui| . 
This is the main conjecture of causal set kinematics to which the previously mentioned length 
conjecture is really just a corollary. Some evidence for the validity of the Hauptvermutung 
has already been established [uj] but a full proof has not yet been formulated. However, 
if the correspondence between a causal set and a manifold is to be (nearly) unique, then 
not only should the length of the longest paths produce the geodesic length, but the longest 
paths themselves should produce the geodesic. 

In this work we look for evidence that the longest paths in faithfully embeddable causal 
sets approach the appropriate geodesic curve as the density of causal set elements increases. 
This work also constitutes an indirect test of the Myrheim length conjecture in curved space- 
times, because we use this conjecture to identify which paths might produce the geodesic 
upon embedding. The details of how this is done are given below. 
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III. COMPUTATIONAL DETAILS 


A. Overall Approach 

If we consider points y t in a spacetime manifold Ai, we can define an interval Im between 
two causally related points y p and yi as the intersection of the future of the pastnrost point 
y p with the past of futuremost point ]jf. Let us also define Sm as the geodesic between y p and 
yf. Then, our general approach is as follows: (a) we generate many faithfully embeddable 
causal sets C = {xi, -<} containing N elements in the form of intervals Ic(xi, xn) by selecting 
points in Im(Uv> 2/f) v i a a Poisson process; (b) in each I c we sample the longest paths Sc 
between X\ and and determine the average goodness-of-ht between the images of the Sc 
in I a/! and Sm] ( c ) we then repeat steps (a) and (b) for larger causal sets and examine the 
behavior of the average goodness-of-fit as N increases. 


B. Random Sprinklings 


As stated previously, we generate faithfully embeddable causal sets by selecting points 
from a manifold via a Poisson process - this is also called a random sprinkling into the 
manifold. Justification for why this should be a Poisson process is given elsewhere 0-0 
The procedure by which we perform these sprinklings is described in 0; but we briefly 
explain it here for completeness. 

The spacetimes used in this study are given by the metric 


ds 2 = Q 2 rj a f)dx a dx 3 , 


( 1 ) 


where fl 2 is the conformal factor and r) a p is the Minkowski tensor. The sprinklings are per¬ 
formed by a (double) rejection method similiar to the methods described in Jf| and Jf|. The 
interval is enclosed in a box and points are randomly selected within this box; if a selected 
point is outside the interval it is rejected, otherwise, it is kept - this is the first rejection. 
In flat spacetime (fl 2 = 1) this first rejection provides the desired uniform distribution of 
points. In curved spacetimes, however, remaining points face a second rejection to ensure 
that the final points are distributed uniformly with respect to the volume form Q d , where 
d is the dimension of the spacetime. Each remaining point is associated with a random 
number w within the range 0 < w < where 0(( iax is the maximum value of the volume 
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form within Im . If w is greater than the value of the volume form evaluated at the point in 
question, the point is rejected; otherwise, it is kept. This process continues until N points 
are sprinkled into /». 

When partially ordered according to their causal relations in Im , these N points, without 
the manifold, comprise a faithfally embeddable causal set interval Ic- At the end of the 
sprinkling process, the selected points are sorted according to increasing value of the time 
coordinate. The corresponding elements in Ic are labeled with a number from 1 to N. 
However, keep in mind that, exclusive of X\ and x n, this natural labeling does not necessarily 
correspond to an ordering of the elements. For example, element X 25 may preceed element 
X 50 (%25 -< X 50 ) or they may not be related in any way. 

C. Finding the Paths 

Having formed the causal set interval, the longest paths between x\ and Xn must be 
found. Starting with element X\ we form the set of elements D 1 linked to its “future”; that 
is, D 1 = {xj G Ic\x 1 + Xj}. For each of these daughters of an we calculate the length L(xj) 
of the path from Xj to Xn using an algorithm due to Sorkin |lOj. We save only those elements 
of D\ for which L(xj ) is maximum. This forms the subset Y D \ of the youngest daughters of 
X\. Every element in Y D \ is in a longest path from x\ to xn- Then, to complete a longest 
path we can choose any element of Yd 1 and repeat the process until we reach an element for 
which the only daughter is Xjv- 

In general, there are many paths that traverse an interval. In fact, the number of paths 
can become quite large. Table 1 lists the average number of longest paths (for 10 trials) 
for various interval sizes N for causal sets generated by sprinkling into 1+1-dimensional 
Minkowski space. As the table shows, the number of longest paths is on the order of 10 10 
for intervals as small as N — 2000. With so many paths in each Ic it is computationally 
impractical to compare the image of every Sc from all the sprinklings to Sm- To circumvent 
this problem we conduct our study using a Monte Carlo style sampling of the paths in each 
trial Ic- 
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Table 1. Average number of longest paths Sc versus interval 
size for ten trial sprinklings. The notation a ( b ) means a x 10 6 . 


N 

(Sc) 

100 

6.70 (1) 

500 

5.18 (3) 

1000 

3.69 (6) 

1500 

1.15 (8) 

2000 

9.97 (9) 


Our procedure is the following: (a) we obtain a trial path ,s'i by selecting the first elements 
(smallest label) in the various Yoy (b) Using the coordinates from which the elements of .Si 
were sprinkled, we calculate the goodness-of-ht of the image of S] to Sm using a y 2 statistic 

1 (2) 


where n is the number of elements in the longest path, d — 1 is the number of spatial 
dimensions of the manifold, the Xj are the spatial coordinates of points in the image of the 
path .s'] in M. and the y are the spatial coordinates of those points along the geodesic having 
the same time coordinate as the point with Xj. So, in standard chi-squared language, the 
points from the path play the role of the “observed” values, the points from the geodesic play 
the role of the “expected” values, and n is the number of “degrees of freedom.” (c) We then 
find a new path by randomly choosing any element, from among all the elements that are 
along any longest path, that has more than one youngest daughter, randomly choosing one 
of the youngest daughters, then completing the path with the first elements of the remaining 
Y Dj . (d) We calculate the y 2 for this new path and save this path if it is a better fit than 
the previous path, otherwise it is discarded, (e) We repeat steps (c) and (d) for a number 
of trials equal to 5 x L(x i). The prescription for the number of trials is adopted so that 
the number of trials automatically scales upward as the number of elements in the causal 
set, and therefore the number of Sj to be sampled, increases; the constant 5 is found to be 
sufficient to produce the best results. The path that emerges from this process is taken as Sc 
for that trial Ic- We conduct multiple trial sprinklings (typically 100 for flat spacetimes and 
1000 for curved) for a given size N and compute the average of the y 2 values for each Sc- 
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In flat spacetimes small enough that all of the longest paths can be checked, this procedure 
sucessfully found the absolute best fit in most of the cases and nearly so in the others jlfij ]. 

It is worth pointing out that there is a sense in which this approach attempts to mimic the 
role of an eventual causal set dynamics. One expects that a sum-over-causal sets dynamics 
should reveal a statistical preference for paths near the geodesic in much the same way 
as the sum-over-histories approach in quantum mechanics reveals a statistical preference 
for classical trajectories. To gauge whether the expected longest path does approach the 
geodesic in the large N limit, we calculate the average chi-squared < y 2 > of all of the 
paths sampled during the Monte Carlo process for all trial sprinklings of a given N. What 
we expect is that, if the conjectures hold, the expected longest path will show consistently 
better fits as N increases. It is important to keep in mind that it is the trend of the fits that 
matters here moreso than the fits themselves. 

In flat spacetimes, it is straightforward to identify the actual geodesic curves Sm- For 
curved spacetimes the Sm are found numerically using the shooting method. The analyt¬ 
ical details of how the geodesic equations are set up for doing these calculations are given 


elsewhere 


M- 


The specific coordinates y 1 - used to calculate the y 2 values in Eq. m 


are 


then found by interpolating between the numerically calculated points. 


IV. RESULTS AND CONCLUSIONS 

A. Flat Spacetimes 

Figure 1 shows the results for 1+1-dimensional Minkowski space. Panel (a) shows a 
uniform sprinkling of 5000 points into an interval together with the geodesic Sm and a best 
longest path Sc determined from the Monte Carlo process described previously. We can see 
from this plot that, modulo statistical fluctuations, Sc produces a reasonably good visual 
fit to Sm- However, it is the decreasing trend of the < y 2 > in panel (b) of this figure that 
clearly establishes the key result that the expected longest path in the interval approaches 
the geodesic as the number of causal set elements increases. Similar results for 2+1- and 
3+1-dimensional Minkowski space are shown in Figure 2 and Figure 3, respectively. For 
those cases (figures 2a, 3b, 3c, and 3d) only projections of uniform sprinklings are shown. In 
those projections the points do not visually appear to be as uniform as in figure la because 
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in the higher dimensional cases there is much more volume near the center of the interval 
than near x\ and xn causing more of the projected points to be centrally located. As with 
the 1+1-dimensional case, we can clearly see in figures 2b and 3a that < y 2 > decreases as 
N increases. These flat spacetime calculations constitute direct numerical evidence in favor 


of the conclusions reached by Brightwell and co-workers j§|. 


B. Curved Spacetimes 

Figure 4 shows the results for a 1+1-dimensional conformally flat spacetime with con¬ 
formal factor Q 2 = ( xt ) 2 . Panel (a) shows a uniform sprinkling of 5000 points. Careful 
inspection of this sprinkling shows that most of the volume (and therefore the points) occur 
at larger values of x and t. Panel (b) shows the key result that the expected longest path 
approaches the geodesic curve as N increases. For a visual sense of what occurs in this 
spacetime, panel (c) compares Sm with Sc also showing the corresponding flat spacetime 
geodesic for comparison. This figure shows that this Sc does appear to follow Sm , especially 
at larger times where there are many more points in the sprinkling. 

Figure 5 shows the results for a 2+1-dimensional conformally flat spacetime with con¬ 
formal factor h2 2 = (x 4 + y 4 ) /t e . Panel (a) shows the projections of a uniform sprinkling 
of 5000 points. These projections show that most of the volume occurs at smaller times 
and larger values of x and y. Panel (b) shows the key result that the expected longest path 
approaches the geodesic curve as N increases. In panels (c) and (d) we see the x — t and y — t 
projections of Sm and an Sc- In a similar way as in the 1+1 case, we see that at smaller 
times, where most of the points are, the best longest path does appear to follow along the 
geodesic curve. 

Figure 6 shows the results for a 3+1-dimensional conformally flat spacetime with confor¬ 
mal factor Q 2 = t 4 e x2 ~ 8y . This conformal factor increases with x and t, decreases with y , and 
has no z dependence. Therefore, the features of the interval are that there is more volume 
at larger values of x and t, smaller values of y, and uniform in 0 . Projections of a uniform 
sprinkling of 5000 points into an interval for this spacetime, shown in panels (b), (c), and 
(d), reveal these features of the interval. Panel (a) confirms that even for this spacetime, the 
expected longest path approaches the geodesic curve as N increases. This fact is supported 
visually by panels (b), (c), and (d) which also show that our best longest path is an excellent 
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fit to the geodesic curve. 


C. Conclusions 

In this paper we have provided computational evidence that the path-geodesic corre¬ 
spondence, implied by key ideas in causal set research, holds for flat and conformally flat 
curved spacetimes. We further suggest that by establishing that the longest maximal chains 
approach the actual geodesic curves in the spacetimes we have considered, we have also pro¬ 
vided evidence in favor of the wider conjecture, the Hauptvermutung, at the core of causal 
set kinematics. 
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V. FIGURE CAPTIONS 

FIG 1. Results for 1+1-dimensional Minkowski space. Panel (a) shows the sprinkled 
points Xj, the geodesic Sm, and one of the best longest paths Sc, determined by the Monte 
Carlo process described in the text. Panel (b) shows the decrease in the < y 2 > versus 
interval size N. The dashed line is just a numerical fit to the data points the purpose of 
which is to emphasize the trend of the data. 

FIG 2. Results for 2+1-dimensional Minkowski space. Panel (a) shows x-t and y- 
t projections of the sprinkled points, the geodesic and its projections, and one of the best 
longest paths and its projections. Panel (b) shows the decrease in the < y 2 > versus interval 
size N. The dashed line emphasizes the trend of the data. 

FIG 3. Results for 3+1-dimensional Minkowski space. Panel (a) shows the decrease in 
the < y 2 > versus interval size N with a dashed line to emphasize the decreasing trend of 
the data. Panels (b) through (d) show projections of the sprinkled points, the geodesic, and 
one of the best longest paths. 

FIG 4. Results for 1+1-dimensional curved spacetime using the metric of Eq. 0 with 
conformal factor U 2 = {xt) 2 . Panel (a) shows the uniform sprinkling of 5000 points into 
an interval. Panel (b) shows the decrease in the < y 2 > versus interval size N. Panel 
(c) compares one of the best longest paths, Sc, to the geodesics, Sm, f° r both the curved 
spacetime and the 1+1-dimensional flat spacetime. 

FIG 5. Results for 2+1-dimensional curved spacetime using the metric of Eq. 0 with 
conformal factor fl 2 = (a; 4 + y A ) / 1 6 . Panel (a) shows x-t, y-t , and x-y projections of a 
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uniform sprinkling of 5000 points. Panel (b) shows the decrease in the < % 2 > versus 
interval size N. Panels (c) and (d) compare projections of one of the best longest paths, Sc, 
to those of the geodesics, Sm, f° r both the curved spacetime and the 2+1-dimensional flat 
spacetime. 

FIG 6. Results for 3+1-dimensional curved spacetime using the metric of Eq. ® with 
conformal factor Q 2 = t 4 e x2 ~ Sy . Panel (a) shows the decrease in the < + 2 > versus interval 
size N. Panels (b) through (d) show projections of the sprinkled points and compares 
projections of one of the best longest paths, Sc, to those of the geodesics, Sm, for both the 
curved spacetime and the 3+1-dimensional flat spacetime. 
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